Gauge Singlet Scalars as Cold Dark Matter 



John McDonald^ 

Cosmology and Astroparticle Physics Group, University of Lancaster, Lancaster LAl 4YB, UK 

February 7, 2008 

Abstract 

In light of recent interest in minimal extensions of the Standard Model and 
gauge singlet scalar cold dark matter, we provide an arXiv preprint of the paper, 
published as Phys.Rev. D50 (1994) 3637, which presented the first detailed 
analysis of gauge singlet scalar cold dark matter. 



j.mcdonald@lancaster.ac.uk 



1 



Abstract 



We consider a very simple extension of the standard model in which one or more 
gauge singlet scalars Si couples to the standard model via an interaction of the form 
AsS^SiH'^H, where H is the standard model Higgs doublet. The thermal relic density of 
S scalars is calculated as a function of As and the S mass mg. The regions of the (ms, As) 
parameter space which can be probed by present and future experiments to detect 
scattering of S dark matter particles from Ge nuclei, and to observe upward- moving 
muons and contained events in neutrino detectors due to high-energy neutrinos from 
annihilations of S dark matter particles in the Sun and Earth, are discussed. Present 
experimental bounds place only very weak constraints on the possibility of thermal 
relic S scalar dark matter. The next generation of cryogenic Ge detectors and of large 
area (lO^m^) neutrino detectors will be able to investigate most of the parameter space 
corresponding to thermal relic S scalar dark matter with ms ^ 50GeV, while a 1 km^ 
detector would in general be able to detect thermal relic S scalars with ms ^ lOOGeV 
as a dark matter candidate and would be able to detect up to ms ^ SOOGeV or more 
if the Higgs boson is hghter than lOOGeV. 

1 Introduction 

There is strong evidence that the mass density of the Universe is mainly composed 
of some non-hadronic form of dark matter^^'^l. Direct observation of galaxies and 
clusters of galaxies'^l indicates that Q = 0.1 to 0.3, where Q is the ratio of the mass 
density to the critical density in the Universe at present. Nucleosynthesis constrains 
the density of hadronic dark matter to satisfyt^' 0.011 < il^h^ < 0.019, where h = 
0.5 — 1 parameterizes the uncertainty in the observed value of the Hubble parameter. 
Inflation and naturalness considerations W suggest that Q — 1. Although it seems 
possible that baryons could just about account for D, — 0.1 dark matter, it would 
not be possible for primordial density perturbations to grow sufficiently in a baryon 
dominated Universe to allow galaxy formationt^l to be consistent with the magnitude of 
temperature fluctuations of the cosmic microwave background radiation as observed 
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by COBE'^I. This requires the addition of a density of non-hadronic dark matter, 
preferably cold dark matter (CDM)'"''!. It would also be difficult to explain, if halo 
dark matter was hadronic in nature, how all the hadrons in galactic halos could be 
hidden'^l. Searches for faint stars support the conclusion that the halo dark matter 
cannot primarily be baryonic'^l (although recent observations of microlensing by dark 
objects in the galactic halo do show that at least some baryonic halo dark matter 
exists'^^). Thus it is likely that the Universe is dominated by a density of CDM 
satisfying I^cdm ^0.1. The age of the Universe imposes an upper limit on fl, flh^ ^ 1 
[^1 . This leaves a window for which a density of particles can consistently serve as the 
primary component of the halo dark matter, 0.025 ^ Jlh^ ^ 1. 

In this paper we will study in some detail an extremely simple and natural extension 
of the SU(3)c X SU(2)l x U(1)y standard model, namely, the addition of one or more 
gauge singlet complex scalars Si. These scalars, if stable, can in principle account for a 
density of CDM. Stabihty of the scalars can most simply be guaranteed if a continuous 
or discrete symmetry exists under which the gauge singlet scalars are the lightest 
particles transforming non-trivially. (Additional continuous and discrete symmetries 
are a common feature of many extensions of the standard model, serving to simplify 
the models and to eliminate phenomenologically unwelcome interactions such as those 
leading to baryon and lepton number violation or to flavor-changing neutral currents.) 
In addition, it is necessary that the Si do not acquire vacuum expectation values, which 
in turn requires that they have positive mass squared terms. This model for CDM is 
essentially determined by just three parameters: the Higgs boson mass mh, the S scalar 
mass mg, and the coupling of the S scalars to the Higgs bosons Xg. In particular, we 
will consider the thermal relic density of S scalars, coming from S scalars freezing-out 
of thermal equilibrium. This is the simplest and most natural origin of a relic density of 
S scalars, although in principle other possibilities exist, such as S scalars originating in 
the out-of-equilibrium decay of some heavy particle. We will be particularly interested 
in the possibility of detecting S scalar cold dark matter, as a function of As, mg and 
mh, either via direct detection of the recoil energy coming from elastic scattering of 
S dark matter particles from Ge nuclei'^^'^^l, or by observing upward moving muons 
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or contained events in neutrino detectors, produced by high-energy neutrinos coming 
from S annihilations in the Sun or in the Earth[i2-i7]_ 

The paper is organized as follows. In section 2 we discuss the thermal relic density 
of gauge singlet scalars in the Universe at present. In section 3 we discuss the elastic 
scattering of S scalars from Ge nuclei. In section 4 we discuss the rate of upward moving 
muons and contained events produced by high energy neutrinos due to S annihilations 
in the core of the Sun and of the Earth. In section 5 we give our conclusions. In 
the Appendix wc give some details of the calculation of the upward-moving muon and 
contained event rates. 



2 S scalar dark matter 

We consider extending the standard model by the addition of terms involving the S 
scalars 

Ls = d^Sld^Si - m^SjSi - AgS/SiHtH (2.1), 

where i = 1, ...,N. This model has a global U(l) symmetry, Si c'"Si, which guaran- 
tees the stability of the Si scalars by eliminating the interaction terms involving odd 
powers of Si and Sj which lead to Si decay. We first consider the case N = 1. In order 
to calculate the relic density arising from S scalars freezing out of thermal equilibrium 
we will use the usual Lee- Weinberg (LW) approximation!^^! to solve the rate equation 
for the density of Si scalars. The rate equation is given by 

^ = -3Hns- < aannVrei > (4 ' 4) (2-2). 

Cann IS the SS^ annihilation cross-section, Vj-ei is the relative velocity of the annihilating 
particles and H is the expansion rate of the Universe. The angular brackets denote 
the thermal average value. (2.2) gives the number density of S scalars ns. The total 
density of S and S^ scalars is then 2ns. The equilibrium S density n^, ioT ms/T ^ 1, 
is given by 



V27rT 



j e-^ (2.3). 
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The approximate solution of (2.2) is then found by rewriting (2.2) as 

df _ < Q-annVrel > ( ^2 f2\ fr, y^^ 

dT~ K 

where f = f|, fo = ^ and K = (47r3g(T)/45M|i)^/l g(T) is the number of degrees 
of freedom with masses smaller than T. In this we are assuming that the number of 
degrees of freedom in thermal equilibrium with the photons, g{T), is constant around 
the S freeze-out temperature Tfs, i.e. no particle thresholds at T Tfs, and also that 
the Universe is radiation dominated. The LW solution is given by assuming that f = f o 
until the temperature at which 



df„ 



dT 



< (JannVrel > ^2 ^2.5) 



K 

is satisfied, which defines the S freeze-out temperature, Tfs. Then for T < Tfs one 
solves (2.4) with fo = on the right hand side and with f(Tfs) = fo(Tfs). The freeze-out 
temperature is then obtained from 

■-l^f ., .T^* ,J (2.6). 



■■'V(1-3xb/2)(2txb)''"/ 

where Xfs = Tfs/ms and A = <£as^vj;ei> ^j^g present total mass density in S scalars 
and antiscalars is then 

^ ^ Ps + Pst _ 2 g(T.) K /Tn (1 - 3xfs/2) 

^ Pe g(Tfs) T^Xfs < (TannVrel > V Pc / (1 " Xfs/2) ^ ' ^' 

where it has been assumed that < CannVrci > is T independent, = 7.5xlO^^'^h^GeV^ 
is the critical closure density of the Universe at present (h = 0.5 to 1) and T^ is the 
present photon temperature. 

In order to calculate < (JannVrei > we need the SS^ annihilation modes. These are 
shown in Figure 1. The corresponding contributions to < (TannVrei > are given by 
SSt ^ h°h° : 



647rmg y mg 



1 - (2-8a), 



sst ^ w+w- 



1 / _ 2m|y \ Ajmt, / _ <\ 

2 1^ m2^y j 87rmi((4mi-m2)^ + m2r2) ^ mi ^ ^ ' ^' 



4 



SSt ^ Z°Z° : 





Xlmi ( ^ ^ ^^"^ 



' ' IGTrmi ((4mi - mgf + m^r^) V " 



SSt ^ ff : 



7rg2 ((4mi-m2f + m2r2) V 
Here the fermion Yukawa coupling is Af = nif/v where v = 250GcV and nif is the 
fermion mass, nih is the Higgs boson mass, and Fh is the Higgs decay width, for which 
we use the standard model values'^^l. We should note that the assumption made in 
arriving at (2.7), that < (JannVrei > is T independent, is strictly true only for freeze-out 
temperatures small compared with the electroweak phase transition temperature Tewj 
where 

2.4mh 

Tew = (2.9), 

+ 0.62 (^"l 



and mt is the t quark mass. The thermal expectation value of the Higgs field is given 

/ \ 1/2 

by < h° >T= V I 1 - ) . Thus for Tfs ^ Tew the effective mass of the W and Z 

V EW / 

bosons goes to zero, while the < h° >t dependent S and h° masses differ from their 
zero temperature values. In practice, however, Tfs ^ Tew occurs only for very large 
S masses ^ ITeV. In this limit, the S mass is essentially determined by the constant 
mass term m in (2.1) and so is effectively < h° > independent, while from (2.8b) and 
(2.8c), in the limit ms ^ mw and ms ^ mh/2 the contribution from S annihilations 
to W and Z bosons reduces to 3 times the contribution from annihilations to the Higgs 
boson (2.8a), as expected in the SU(2)xU(l) symmetric limit. This is T independent. 
A large Higgs boson mass doesn't alter this conclusion, since from (2.9) Tew is of order 
the Higgs mass, and so if Tfs is of order Tew, then ms, which, as discussed below, is 
much larger than Tfs, will also be much larger than mh. As a result, we can neglect 
mh in the propagators of (2.8). Thus in practice we can use the cross-sections (2.8) 
calculated with the T=0 value for < h° >, < h° >= 250GeV. 

Using these contributions to < CTannVrei > we solve (2.6) self-consistently for the 
freeze-out temperature and then obtain from (2.7) the resulting dark matter density. 



In Figures 2(a)-2(d) we give plots of the dark matter density as a function of ms and As 
for various values of nih- In these we show the contours for Qsh^ = 1.0, corresponding 
to the upper limit from the age of the Universe, fish^ = 0.25, which is the smallest 
value for which a critical density (f^s = 1) of S scalars can occur, and Osh^ = 0.025, 
corresponding to the smallest value for which S dark matter could make up the primary 
component of the galactic halo. In Table 1 we give values of the freeze-out temperature 
for various values of As, ms and nih. Typically ms = (10 — 30)Tfs for the range of 
parameters we are considering. We have assumed mt = 120GeV throughout. We find 
that increasing the t quark mass to 200GeV makes only a very small difference to the 
results. 

Table 1. S freeze-out temperature, x^^ = ms/Tfs 



As 


ms 


mh 




\s 


ms 


mh 


XfS 


1 


30 


60 


34.5 


0.01 


30 


60 


25.5 


1 


30 


100 


23.5 


0.01 


30 


100 


14.6 


1 


30 


300 


18.5 


0.01 


30 


300 


9.6 


1 


30 


500 


16.5 


0.01 


30 


500 


7.7 


1 


100 


60 


27.9 


0.01 


100 


60 


18.9 


1 


100 


100 


28.3 


0.01 


100 


100 


19.3 


1 


100 


300 


26.5 


0.01 


100 


300 


17.5 


1 


100 


500 


23.9 


0.01 


100 


500 


14.9 


1 


1000 


60 


25.1 


0.01 


1000 


60 


16.1 


1 


1000 


100 


25.1 


0.01 


1000 


100 


16.1 


1 


1000 


300 


25.1 


0.01 


1000 


300 


16.1 


1 


1000 


500 


25.2 


0.01 


1000 


500 


16.2 



Prom Figures 2(a)-2(d) we see that for large values of As (larger than 0.1), which 
is particularly interesting from the point of view of the phenomenology of S scalars, 
in order to have a density of S scalars which can account for a critical density of dark 
matter (Osh^ ^ 0.25), we require the S mass typically to be ^ lOOGeV or ^ 500GeV. 
More generally, as the Higgs boson mass increases, the value of As for which S particles 
of mass less than about lOOGeV can account for dark matter increases, from about 
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As = 0.01 - 0.1 for nih = 60GeV to As > 1 for nih > 300GeV. This could be important 
with respect to the possibihty of producing S particles via Higgs decay at future multi- 
TeV hadron colliders such as the CERN Large Hadron Collider (LHC) '^^1. 

These results are for the case N=l. If we consider several scalars of equal mass and 
couphng strength to the Higgs (for example, if the Si were a multiplet under a global 
symmetry or even gauge symmetry with a symmetry breaking scale large compared 
with mw), then it is easy to see that the total density in Si and Sj is just the sum 
over each individual Si density, since each Si annihilates only with its own antiparticle. 
Thus 

QS Total = Si%=NQs (2.10) 

Since Qs cx - — - — - oc t^, we see that Figure 2 still holds if we replace As by 
As — As/-\/N on the horizontal axis. Thus for a given value of ilg and ms the value of 
As is increased by -\/N. This will increase strength of interaction with matter and so 
the observabihty of S dark matter for N > 1. 

3 Elastic scattering of S dark matter particles from 
nuclei and constraints from Ge detectors. 

In this section we consider the constraints on As and ms following from direct detec- 
tion of S dark matter particles via elastic scattering of S scalars from Ge nuclei^^^'^^l. 
It will be assumed throughout that S dark matter accounts for the halo dark matter 
density. Although the simplest possibihty for the origin of a relic density of S par- 
ticles is from S freeze-out, in principle there are other possibilities. For example, if 
a heavy particle such as a heavy right-handed neutrino N decays to S particles (via 
the Higgs-mediated process N — > i^lS^S in the case of right-handed neutrinos) at a 
temperature below the S freeze-out temperature (typically between IGeV and 50GeV 
for 20GeV ^ ms ^ lOOOGeV) then the S particles so produced will not return to an 
equilibrium density and will result in a relic S density different from the thermal relic 
density. In this case halo dark matter could, in principle, be accounted for by any com- 
bination of ms and As. Thus it is important to consider generally what constraints on 
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the parameters of the model are imposed by experimental observations, as well as to 
compare the constraints with the thermal relic density as a particular example. 

The S scattering cross-section from quarks via Higgs exchange gives an effective 
interaction 

^eff = ^StSqq (3.1). 



m-^ 

Using the expression for the nuclear matrix element^^^'^^' 

2 — 

< N|Eq mqqq|N >= (7.62)— mN^'N^'N , 
we see that the effective interaction with a nucleus is given by 



Leff = (7.62)|-^StSV'NV'N (3.2) 



and that the cross-section for coherent S-nucleus scattering is given by 

-s-N = (7.62)^ .J... A (3.3). 

(277r) (ms-FmNj 

In general cts-n must be multiplied by a correction factor (^N(nis), which accounts 
for the fact that at large enough momentum transfer the scattering ceases to be a 
coherent scattering with the whole nucleust^"^'. Wc will use a correction factor based 
on integrating a Gaussian nuclear form factor over the Maxwellian velocity distribution 
of the halo dark matter particles '^^'^^^ : 



0.573 

CN(nis) = ^— 



cxp(-(Yjb))erf(^^ 
;i + b) erf(l) 



(3.4), 



where 



b = gV Wge(^^^^^)2 ' (3-5) 

rcharge = 5.1(0.3 + O.SQA^/^^GeV-^ , 

and V = V3oo300kms~^ is the mean halo velocity dispersion of the S particles, which 
is related to the galactic rotation velocity in the isothermal sphere model, Vroti by 



V — \ j^l2vrot- For the case of scattering from Ge we find that the full cross-section is 
given by 

_„^„_36 _2 CN(ms) ^ lOOGeV y^AsV 
(7s_Ge = 5.7x10 cm -2 ?rr (3.6 , 



76GeV ) 



where CN(nis) is given by (3.4) with 

b = 2.2v^oo/(l + 76GeV/ms)' 



In order to compare with experiments we need the rate of interaction of the halo dark 
matter particles with a detector per kg per day. This is given by (without energy 
threshold) [23] 

P / 8 rj,yp,,as-NvA^s) . n /100GeV\ /100GeV\ , 

R = U~ ^ 0.069 po.4 V300 0-36 kg d 

KSnJ nismN V mN / V ms / 

(3.7), 

where po.4 is the density of S scalars in the halo (ph) in units of 0.4GcV cm"^, ?7v ~ 1.3 
is a correction for the motion of the Sun and the Earth, and cr^ is the corrected S-N 
cross-section in units of 10~^^cm^. 

It should be noted that the correction factor (^N(iris) is not accurate for dark matter 
particle masses much larger than lOOGeV t^^l. However, we will see that the exper- 
imental constraints in the present model are most important for S masses less than 
about lOOGeV, in which case the correction factor (3.4) is accurate to about 10% ^^1. 

In Figure 3 we show the event rate as a function of ms and As for the cases 
mh = 60GeV, lOOGeV and 300GeV. We also show the contours corresponding to the 
thermal relic S dark matter region of the parameter space, 0.025 ^ f^s ~ 1- (We have 
assumed po.4 = V300 = 1 throughout). 

The present experimental upper bound on R corresponds approximately to 100kg~^d~^ 
for ms ^ lOGeV [^^'^^'^^l. In general, present ionization detectors may be able to achieve 
a sensitivity of about 10kg~^d~^ l^^l, while in the future cryogenic Ge detectors (such 
as a proposed SOOg^'^Gc + SOOg^^Gc detector'^^1) should be able to achieve a sensi- 
tivity of 0.1kg~^d~^. We see from Figure 3 that in order to constrain the thermal 
relic S region of parameter space we need an upper bound on R which is less than 
100kg~^d~^. For an upper bound on the cross-section of 10kg~^d~^, we can probe 
a small region of the thermal rehc parameter space corresponding to As ^ 0.06 and 
ms ^ 20GcV. In order to significantly constrain the possibility of a critical density of S 
dark matter, 0.25 ^ f^sh^ ^ 1-0, we require R ^ lkg~^d~^. R ^ 0.1kg~^d~^ would allow 
us to detect or exclude almost all thermal relic S dark matter for ms ^ 50GeV, while 
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R "Z 0.01kg~^d~^ would detect almost all thermal relic S dark matter for ms ^ lOOGeV. 
These conclusions for ms ^ lOOGeV are essentially independent of mh, as can be seen 
by comparing Figures 3(a), 3(b) and 3(c). For ms ^ lOOGcV, the amount of thermal 
relic parameter space which can be experimentally searched decreases as m^ increases. 
From Figure 3(a) we see that even with mh as small as 60GeV, in order to constrain 
the thermal rehc dark matter for ms ^ lOOGeV we would need R ^ 0.01kg~^d~^, and 
even if an upper bound as low as O.Olkg^^d^^ could be achieved, this would not be 
sufficient to constrain the possibility of a critical density of thermal relic S dark matter 
for ms > lOOGeV. 

Thus we can conclude that since present Ge ionization detectors give an upper 
bound on R of about 100kg~^d~^ (and are not expected achieve a sensitivity better 
thatlOkg^^d^-*^), present attempts at direct detection of dark matter can at best only 
impose a very weak constraint on the possibility of thermal relic S dark matter. The 
next generation of cryogenic detectors should be able to effectively search for S dark 
matter up to at least 50GeV. Thermal rehc S scalars significantly heavier than lOOGeV 
are probably beyond the reach of future Ge detectors, even if the Higgs mass is as small 
as 60GeV. 

4 High energy neutrinos from SS^ annihilation in 
the Earth and the Sun. 

In this section we calculate the fiux of upward-moving muons and the rate of contained 
events in neutrino detectors due to high-energy neutrinos (> 2GeV) resulting from 
annihilations of SS^ pairs in the core of the Earth and the Sun'^^"-*^^! . 

The rate of upward-moving muons at the surface of the Earth due to annihilations 
in the Sun is given by '^^^ 

Tdetector = 1.27xlO-29Cm2SiaibiSFBF < Nz2 >F i m-^yr-' (4.1), 

where C is the capture rate in the Sun in units s~^, ai and bi are the neutrino-scattering 
and muon-range coefficients, summed over i = i/^ and (a^,^ = 6.8, aj?^ = 3.1, 
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"dN" 


/ 


dE 



EME (4.2), 

Fi 



^. is the differential energy spectrum of neutrino i at the surface of the Sun 



h^^ — 0.51, bx7^ — 0.67) and the Bp are the branching ratios for SS^ annihilations to 

gauge boson, Higgs boson and quark pairs. < Nz^ >pi arc the second-moments of the 
spectrum of neutrino type i from final state F scaled by the S mass squared, 

< Nz^ >Fi= — 

m; 

where 

or Earth resulting from injection of particles in final state F at the centre of the Sun 

or Earth. For the case of annihilations in the Earth one multiplies (4.1) by 5.6x10^, 
corresponding to the ratio of the distance squared to the Sun to the radius squared of 
the Earth I^^]. 

The capture rate is given hy^^'^'^^'^^''^'^^ 

C = C— ^^^SN(T4oFN(ms)fN0NSN/niN (4.3), 

nisvsoo 

where is the S- nucleus elastic scattering cross-section in units of lO^^'^cm^, c = 
5.8xl0^'^s^-'^ for the Sun and c = 5.7xl0"^^s~^ for the Earth. The sum is over all species 
of nuclei N in the Earth or Sun. 0n and fN are given in Table A.l of ref.l7. Sn is 
a factor which takes into account the fact that the S dark matter particle must lose 
sufficient momentum to be captured. For Sn we have 



V2 

esc 

V2 



0i (4.4), 



2 (ms - mN)^ 

where Vesc is the escape velocity for the Sun or Earth (618kms"~^ and llkms"^ respec- 
tively). This has the correct behaviour for An large and small compared to 1 I^^'^^l. 
FN(mg) is a factor which takes into account form factor suppression. The branching 
ratios Bp, corresponding to the rates for S annihilation in the limit of zero relative 
velocity, are directly obtained from (2.8). 

The rate (4.1) assumes that the accretion of S particles by the Sun or Earth and 
their subsequent annihilation are in equilibrium, in which case the annihilation rate 
is given by Fann = C/2. The condition for this to be true is that the age of the solar 
system t© be large compared with the time for equilibrium to be established ta, which 
is defined below. In general C in (4.1) should be replaced by^^^'^^l 

C^Ctanh2(^) (4.5), 
Ta 
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where 

TA = 1/(CCa)'/' ; CA=<av>V2/V2 (4.6). 

< (TV > is the spin-averaged total annihilation cross-section times the relative velocity 
in the limit of zero relative velocity, which can be obtained from (2.8). The effective 
volumes Vj are given by 

Vj = 6.5xlO^^(jms/10GeV)-^/^ cm=^ (4.7a), 

for the Sun l^^l and 

Vj = 2.0xl02^(jms/10GeV)-^/2 cm=^ (4.76) 

for the Earth t^^l. 

A second assumption in obtaining (4.1) is that the capture rate is primarily due to 
single collisions with nuclei ("optically thin" hmit). However, for the case of capture 
due to scattering from iron in the Earth, it has been pointed out that multiple collisions 
can enhance the capture rate '^^1 . The enhancement factor is given by 

exp (Teff - 1) , . 

a{T) = , Teff = T/3_ , (4.8) 

Teff 

where 

T (Js-Fe/ (2.3 X 10"^^ Cm^) 

is the optical depth of the Earth, = 4msmFe/ (mg — mpe)^, and Tg// is the effective 
optical depth of the Earth taking into account multiple collisions. This expression 
is valid so long as Max(l, ln/3_) ^ 6/Tes and /3_ ^ 20; otherwise the enhancement 
must be evaluated numerically, although the largest enhancement occurs typically for 
P_ ^ 20 [^■''1. Wc have included a(r) from (4.8) in our calculations over the range 
where it is valid (and where enhancement is expected to be most important), in order 
to indicate the importance or otherwise of multiple collisions. In practice, we find that 
no enhancement of the event rate in detectors occurs over the range of parameters we 
are considering. 

The < Nz^ >Fi arc related to the muon neutrino and anti-neutrino energy spectra 
coming from annihilation of S in the Sun and Earth, including in the case of the Sun 
the effects of the interactions of the annihilation products and neutrinos with the solar 
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mediuml^^l . For ms > mw the dominant contributions to the < Nz^ >Fi are from 
annihilation to gauge boson and t quark final states, while for ms < m-w the dominant 
final states contributing to < Nz^ >Fi are b quark pairs and possibly Higgs boson pairs 
if mh < mw- In the Appendix we discuss the values of < Nz^ >Fi coming from the 
different final states. 

In order to calculate the capture rate for the case of the Earth one can simply use 
(4.3), since in this case the form factor suppression is small for most values of ms 
and so we can take FN(ms) = 1 '^^1. [Capture is dominated by low momentum transfer 
scattering excpet for ms close to the mass of the scattering nucleus, in which case the 
form factor suppression can be more significant (a factor of 0.72 for the case where 
ms = mpe '^^')-] For the case of capture by the Sun, form factor suppression cannot 
be neglected, making the calculation of the capture rate more complicated. A simple 
expression for the capture rate in this case has been given by Kamionkowski'^^l , which 
is accurate to 5% for dark matter particle masses greater than a few GeV and less 
than a few TeV (see also Ref.[14]). In terms of us-n this may be written as 

TT (ms + mN)^ f ( \ (Aa\ 

C = (7S-N fs(nis) (4.9), 

4 mgm^ 

where 

fs{ms) = 2.04 X 1038exp[-0.0172(m<j - 10)], ms < 80 GeV 

6.10 X 1037(^^/80)-106-0.38[(ms-80)/920]V2^ gQ <ms< 1000 GcV 

1.72 X 1036(m5/1000)-l•«^ ms > 1000 GeV 

(4.10). 

Here C is in s~^ and all masses are in GeV. 

In Figures 4(a)-4(c) we show the results for S annihilations in the Sun for the cases 
mh = 60GeV, lOOGeV and 300GeV respectively, and in Figures 5(a)-5(c) we show the 
corresponding results for the case of annihilations in the Earth. Comparing Figures 
4 and 5, we see that for an upper bound on Fdetector corresponding to the 1MB upper 
boundPs], 

Tdetector < 2.65x10 ^m ^yr ^ (curve (a) in Figures 4 and 5), the strongest 
constraints on the parameter space come from Earth S annihilations, while for upper 
bounds less than around 10~^m^yr~^ the solar S annihilations become more important. 
Thus we will compare the thermal relic parameter space with the Earth S annihilation 
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constraints for the case of the present 1MB upper bound and with the solar annihilation 
S constraints for the case of the bounds expected from future neutrino detectors. 

From Figure 5(a), corresponding to Earth S annihilations with mh = 60GeV, we 
see that at present the 1MB constraints can exclude only a small region of the thermal 
relic parameter space (corresponding to the iron 'resonance' at ms ~ 56GeV). Prom 
Figure 4(a), we see that for an upper bound on Fdetector of lO^^m^^yr^^, upward- 
moving muons from solar S annihilations can exclude a small region of the thermal 
relic parameter space corresponding to As ^ 0.1 and ms ^ 50GcV. [Wc see from 
Figure 5(a) that for this case the bounds due to neutrinos from the Earth are stronger 
that those due to solar neutrinos for ms between 20 and 80 GeV, and can probe a 
significant region of the thermal relic parameter space for ms between 50 and 70 GeV.] 
With an upper bound Fdetector < 10~^m~^yr~^, most of the thermal relic parameter 
space in Figure 4(a) corresponding to fls^^ ^ 0.25 for ms ^ 400GcV and a significant 
region of the thermal relic parameter space corresponding to a critical S denisty for 
ms "Z 50GeV can be investigated, while an upper bound Fdetector < lO^^m^^yr"^ would 
probe the whole thermal relic parameter space up to ms ~ 500GeV. For larger mh, 
the conclusions for ms ^ lOOGeV are essentially unchanged, while for ms ^ lOOGeV 
the amount of thermal relic parameter space which can be investigated for a given 
upper bound on Fdetector decreases as mh increases. 

In order to see how the 1MB upper bound could be improved in the future, we 
can make a rough estimate of the bound which could be imposed by building neutrino 
detectors of larger area. The 1MB bounds follow from a detector area of 400m^ and 
exposure of about 1 year, corresponding to an upper bound of less than about 10 
upward moving muons per year. Following ref . [27] , we can estimate the area of detector 
required in order to achieve a given sensitivity by the area needed to detect one upward- 
moving muon event per year. At present, the MAGRO detector at Gran Sasso, with 
an area ~ lO^m^, is beginning operation ^^^1, while a number of detectors with an 
effective area of order lO'^m^ arc under development (DUMAND^^l; AMANDA l^o]; 
NESTORt^^l). In addition, it has been suggested^^l that a Ikm^ detector is needed 
to observe muons from neutralino dark matter in the GeV-TeV mass range, and that 
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this could be constructed at a cost of order 100 million U.S. dollars '^^1. We see from 
Figure 4 that a detector of area 10"^m^ should be able to probe the region of parameter 
space corresponding to curve (c), which will rule out much of the parameter space 
corresponding to thermal relic S dark matter with ms ^ 50GeV. This conclusion is 
essentially independent of mh, as can be seen by comparing Figures 4(a), 4(b) and 
4(c). For the case of a Ikm^ detector, the area of parameter space under curve (e) 
in Figure 4 could in principle be searched. This would probe the entire thermal relic 
dark matter region for ms "Z l.STeV (SOOGeV, lOOGeV) for the case of mh = 60GeV 
(lOOGeV, 300GeV). Thus in general the thermal rehc dark matter parameter space 
can be probed for ms at least up to lOOGeV. 

We therefore conclude that at present the 1MB upper bound on the flux of upward- 
moving muons can impose only a slight constraint on the possibility of thermal relic S 
dark matter, while many of the S dark matter possibilities with ms "Z 50GeV should be 
within the reach of experiments with an area O(10^)m^ in the not-too-distant future. 
In the more distant future large (Ikm^) detectors should be able detect or exclude S 
dark matter for ms up to at least lOOGeV. 

We also note that comparing the next generation of Ge detectors, which might 
reach scattering rates 0.1kg~^d~^ (Figure 3, curve e), with the next generation of neu- 
trino detectors, which might reach upward-moving muon rates 10~^m~^yr~^ (Figure 
4, curve c), we see that for ms ^ 80GeV the constraint from upward- moving muons 
are dominant, while for ms ^ 80GeV Ge detectors impose a stronger constraint. 

Up to now we have only considered the rate of upward-moving muons in discussing 
constraints on the (As,ms) parameter space. The upward-moving muon flux is the 
most important signal for high-energy neutrinos from the point of view of future large- 
area neutrino detectors, which are specifically designed to detect this flux. However, 
in discussing the present bounds on the flux of high-energy neutrinos due to SS^ 
annihilations, we should also consider the possibility that a high-energy electron or 
muon neutrino could undergo a charged current interaction within the volume of the 
detector'^^l ('contained event'). For the case of neutrinos from the Sun the rate of 
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contained events per kiloton due to electron and muon neutrinos is given by'^^l 

rdetector = 3.3xlO-'^CmsEiaiEpBF < Nz >F i kt-Vr-' (4.11), 

where i is summed over the electron and muon neutrino and anti-neutrino. In the 
Appendix we discuss the values of < Nz >f i coming from the various final states. In 
Figure 6 we show the results for S annihilations in the Sun for the cases mh = 60GeV, 
lOOGeV and 300GeV, while in Figure 7 we show the corresponding results for the 
case of S annihilations in the Earth. The present upper bound on the rate of electron 
and muon contained events in the Frcjus detector is'^'''^^! rdetector < 6.4 kt~"^yr~-'-, 
corresponding to curve (a) in Figures 6 and 7. Comparing with the upward-moving 
muon bounds from Figures 4 and 5, we see that at present the contained event rate 
imposes constraints on the parameter space which are in general weaker than those 
coming from the upward- moving muon flux, although at small ms, mg ^ 20GeV, the 
constraints become similar (and slightly stronger for the case of solar neutrinos). 

So far in this section and in the previous section, we have considered the case of 
just one S scalar. For the case of N scalars of equal mass and coupling, the density 
of each scalar Si contributes a proportion 1/N of the total halo density. The capture 
rate of S dark matter in the Sun and the Earth and the rate of elastic scattering in Ge 
detectors are proportional to the density of Si in the halo times the cross-section for 
scattering from nuclei in the Sun and Earth. Thus the contribution to the event rate 
in a detector is reduced by a factor 1/N for a given Si. The total rate from summing 
over i for a given ms and As is therefore unchanged. However, for a given value of 
the thermal relic density, the value of As for a given ms is increased by a factor V^, 
leading to an increase in the scattering cross-section and so to an increase in the event 
rate in Ge detectors and in neutrino detectors by a factor N for a given thermal relic 
density, thus making the dark matter easier to detect. 

5 Conclusions 

The extension of the standard model by the addition of a gauge singlet scalar provides 
a canonically minimal extension of the standard model which can potentially account 
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for dark matter. It is important, therefore, to consider in some detail the question of 
the rehc density of the gauge singlet scalars and their possible observable signatures. 
In general, present experiments based on observing elastic scattering of halo dark 
matter patricles from Ge nuclei, and on observing upward-moving muons at the Earth's 
surface, coming from muon neutrinos due to dark matter particle annihilation in the 
Sun or the Earth, can only place very weak constraints on thermal relic S dark matter, 
and cannot constrain the possibility that thermal relic S dark matter could account for 
a critical density of dark matter (fig = 1). However, the next generation of cryogenic 
Ge detectors (which hopefully should achieve bounds on the Ge scattering rate of 
0.1kg""^d~^) and neutrino detectors (with an effective area lO^m^) will be able to 
investigate most of the parameter space for thermal relic S scalar dark matter with 
ms ^ 50GeV, while a Ikm^ neutrino detector, as suggested in order to search for 
heavy neutralino dark matter, would be able to exclude thermal relic S dark matter 
for ms ^ lOOGeV (as would a cryogenic Ge detector if it could achieve a sensitivity of 
O.Olkg-M-^) . For a light Higgs mass, equal to 60GeV (lOOGeV) , a Ikm^ detector could 
also exclude heavier thermal rehc S dark matter up to l.STeV (SOOGeV). In general, 
the next generation of cryogenic detectors will be the most effective in searching for 
S dark matter with ms ^ 80 GeV, while for ms ^ 80 GeV the next generation of 
neutrino detectors will be most effective. 

The couphng of a gauge singlet scalar to the standard model Higgs doublet is unique 
in form and inevitably will be a feature of many particle physics models beyond the 
standard model. We believe the results presented here may generally be useful in the 
study of such models and of their cosmological consequences. 

Appendix. < Nz^ >Fi and < Nz >Fi from SS^ annihi- 
lations. 

In this Appendix we give the dominant contributions to < Nz^ >Fi and < Nz >Fi for 
the gauge boson, Higgs boson and quark final states coming from SS'^ annihilation. We 
will use the discussion of Ritz and Seckel f^^l (RS) for the case of the quark final states. 
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while for the case of the gauge boson final states we will follow Ref.[16] and consider 
< Nz^ >Fi and < Nz >Fi to mostly originate from the highest-energy "semiprompt" 
W and Z decays to neutrinos. For the case of the Higgs boson final states we will 
adapt the results of RS to obtain < Nz^ >Fi and < Nz >Fi. 

SSt ^ WW, zz 

In this case the dominant contribution to < Nz^ >f comes from muon neutrinos 
originating in the decays W"*" I-I'^t^ix and Z° ^^Pfx- The mean energy squared of 

2 

the neutrinos is Eq = ^(1 + /5^/3), where (3 is the velocity of the decaying W or Zt-*^^! 
[(3 — (1 — m^/mg)^/^ for the case of the W]. This assumes that in the rest frame of 
the W, the W decays isotropically to final states each of energy mw/2. The branching 
ratio of W"'' to decays is given by 1 divided by the number of SU(2) doublets to 
which W can decay, which gives 1/9 for W decaying to all lepton doublets and 1st 
and 2nd generation quark doublets. Thus, noting that N is the number of neutrinos 
produced per injected boson or fermion pair l^^l, we see that the < Nz^ >Fi following 
from annihilation to W pairs can be estimated to be 

< Nz^ >w ^ ^(1 + /^VS) = 0.028(1 + (3^/3) (Al). 

For the case of annihilation to Z pairs, the branching ratio for Z — > V^Vn is 0.066 '^^1, 
and so the < Nz^ >Fi can be estimated to be 

< Nz' >z 2(0.066)^(1 + /5V3) = 0.033(1 + (3^ /?>) (A2), 

where the factor 2 occurs because either of the Z's produced by S annihilation can 
lead to a z/^. The same results are obtained for \ — V^. For the case of the < Nz >Fi 
one obtains in the same way, for i = e and /x, 

< Nz >w u,^W^ 0.056 (A.3) 
9 2 

and 

< Nz >z^i« 2(0.066)i = 0.066 (A4), 
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where we have replaced the mean squared energy of the neutrinos ^(1 + /3^/3) in 
(A.l) and (A. 2) by the mean squared energy in the rest frame ms/2. The values of 
< Nz >Fi for Vi are equal to those of Vi. 

These results are true for the case where interactions of the W, Z and neutrinos 
with the Sun or Earth are ignored. This is justified for the Earth, but for the case 
of the Sun there is an additional suppression factor due to the absorption of muon 
neutrinos (due to charged current interactions) and loss of neutrino energy (due to 
neutral current interactions) as the neutrinos pass through the Sun t^^l . (The W and Z 
will decay fast enough that the effect of their interaction with the solar medium prior 
to their decay can be ignored '^^1 .) The suppression factors are given by '^^1 

Pi = l/(l + Eori)"+"' (A5), 

where Eq is the initial neutrino energy and n = 2(1) for the case of < Nz^ >Fi (< 
Nz >Fi). a^^ = 5.1, av^ = 9.0, r^,^ = LOlxlQ-^GeV-^ and r^^ = S.SxlQ-^GeV-^ for 
i = e, /X. The unsuppressed < Nz^ >Fi and < Nz >Fi are multiplied by the Pi in order 
to obtain the true < Nz^ >Fi and < Nz >Fi for the case of neutrinos from the Sun. 

It is important to note that the assumption that the < Nz^ >Fi and < Nz >Fi are 
dominated by the "semiprompt" decays of the W and Z is well justified for the case 
of the unsupressed < Nz^ >Fi and < Nz >Fi which is appropriate for the case of 
neutrinos from the Earth. However, for the case of neutrinos from the Sun, because 
the the higher-energy neutrinos from semiprompt decays arc preferentially absorbed 
relative to the lower-energy neutrinos coming from secondary decays t^^l (such as W's 
decaying to pairs of quarks which subsequently decay to neutrinos), the secondary 
decay neutrinos may become important at large S masses. At the end of this Appendix 
we make an estimate of the importance of the secondary decays for the case of the 
Z boson final state, where it is shown that the primary decays dominate < Nz^ >-pu 
(< Nz2 >Fp) for ms up to at least 1.4 TeV (2.2 TeV), and up to at least 860 GeV 
(1.3 TeV) for < Nz >-pi, (< Nz >ff)- Prom the Pigures we see that an underestimate 
of the < Nz^ >Fi or < Nz >Fi by a factor of 2 will make very little difference to 
our conclusions. Thus we expect that in general our results for the case of solar S 
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annihilations will be reliable for ms up to at least ~ 1.5 TeV for the upward- moving 
muons and up to at least ~ 1 TeV for the contained events. 

SS^ ^ tt , bb 

In the case of quark final states, one must consider the details of hadronization and 
fragmentation of the final state quarks, which will produce hadron jets. RS I^^l have 
used the results of the Lund Monte Carlo program, which simulates the final states of 
e"*", e~ annihilations into fermion pairs, in order to calculate the values of < Nz^ >Fi 
and < Nz >Fi due to dark matter particles annihilating to fermion pairs. For the case 
of non-interacting final state quarks (appropriate for S annihilations in the Earth) , one 
can use the RS results directly. In general, the < Nz^ >f i and < Nz >f i are given by 

[13] 

< Nz2 >F i= I < y2 > (< z| > -iz^) {A.6) 

and 

N 

< Nz >Fi= - < y >< zf > (A7), 

where z|[ = m^/mg, N, < y° >, the hadron mass mn, and < Zp > are given in Tables 
2 and 3 of Ref. [13]. In this we have assumed that the mean hadron energy (scaled 
by the S mass) whrn the hadron decays, < zh >, is equal to the hadron energy after 
fragmentation < zf >, which is true if the hadrons are not slowed by the astrophysical 
medium (Sun or Earth) before they decay. We will show below that this is in general 
justified for the case of interest to us here. One has to correct (A6) and (A7) for the case 
of ms near the thershold for producing a hadron, since in this case energy conservation 
implies < zp >— * 1. RS make the replacement < Zp >^< Zp > +(1— < Zp >)z^ 
in order to take this into account ^^^1. Using (A6) and (A7) (corrected for thresholds) 
and the results of Ref. [13] we obtain, for the tt final state, 

< Nz2 >t,^= 1.7 X 10-2(1 - 0.04z2,) (A.8) 

and 

< Nz >t,,= 4.7 X 10-2(1 + 0.41z^) (^.9), 
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where i — e or /j,. The same results are obtained for the antineutrinos. For the bb final 

state we obtain 

< Nz^ >b^^= 6.5 X 10-2(1 + 0.39z^) (^.10) 

and 

< Nz >b^.= 2.8 X 10-2(1 + 0.41z^) {A.ll). 

These results are for the case where the interactions with the astrophysical medium 
are ignored. For the case of solar annihilations one has to consider the possible effects 
of hadrons slowing before they decay, as well as the effect of neutrinos losing energy 
or being absorbed as they pass through the Sun. In fact, we can ignore the effect of 
hadrons slowing for the case of interest to us here. For the b quark final state, the 
effect of slowing is only important for Eb > = 470GeV t^^' . But the b quark final 
state is important only when the W and Z final states are kinematically disallowed, 
ms < mw, in which case we can ignore the slowing of the hadrons. For the case of 
the t quarks, one has E^ = (mt/mb)^/2Eb = 2.3TeV for mt = 120GeV. For ms large 
compared with mw, the branching ration to the W final state is much larger than that 
to the t quark final state. [From (2.8) we find Bww/B^^. w 2ms/3mj in the limit of 
large ms.] Thus we see that for values of ms for which the slowing of t quarks becomes 
important (greater than 1 TeV), we can ignore the t quark final states. Therefore in 
general we can ignore the effect of quarks slowing before they decay. 

In order to take account of the interaction of the neutrinos with the Sun, we use 
the method of RS. We simply integrate the differential energy spectrum, including the 
Pi factors from (A. 5): 



<NZ° >FiA= / 



dN 
d^ 



Z (IZ 



Pi (1 + z/zsi) 

where zgi = 1/Tims and < Nz° >FiA is the moment of the neutrino distribution 
including the effect of interactions with the Sun. For {n + ai)z/ zsi small compared 
with 1, the denominator can be expanded to give 

< z°+^ > n + tti 
< z° > Zsi 

Using Table 3 and Eq.(32) of Ref.[13] we find that (< z >, < z^ >, < z^ >) equals 



< Nz"" >FiA=< Nz° >Fi 



(A13). 



(0.13, 4.4 X 10"^ 2.1 X 10-2) for the t quark final state and (0.13, 2.9 x 10"^, 9.5 x 10"=^) 
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for the b quark final state. Thus we obtain, for the case of the t quark final state, 

< Nz^ >t,^ A=< >t,^ [1 - ms/(290 GeV)] (^.14), 

< Nz^ >tj,^ A=< Nz2 >tp^ [1 - ms/(492 GeV)] (^.15), 

< Nz >t,, A=< Nz >t,, [1 - ms/(478 GeV)] (A.16), 

< Nz >tv, A=< Nz >tv, [1 - ms/(765 GeV)] (^.17), 
and, for the b final state, 

< Nz^ >b., A=< Nz^ >b,^ [1 - ms/(422 GeV)] (A.18), 

< Nz^ >bp^ A=< Nz2 >bp^ [1 - ms/(716 GeV)] (A.19), 

< Nz >b.^ A=< Nz >b., [1 - mg/(740 GeV)] (A.20), 

< Nz >bp^ A=< Nz >bp, [1 - ms/(1200 GeV)] (A.21). 

These should be accurate so long as the suppression factors are not too small com- 
pared with 1. However, for the case of the t quark we see that for < Nz^ >tt'^ a the 
approximation breaks down for ms larger than about 250 GeV. In this case an alter- 
native method for estimating the suppression of neutrinos must be used. Prom Table 
3 of Ref.[13] we see that the effect of fragmentation of the t quark is quite small, with 
< zp >= 0.87 and < Zp >= 0.78, compared with 1 for the case without fragmentation. 
In addition, most of the neutrinos come from the primary decay mode to neutrinos, 
t —>■ hfx'^i'n l^^l. This can be seen by comparing the naive estimate based on this decay 
mode with the results of (A. 8) and (A. 9). Assuming that in the rest frame of the 
decaying quark the decay is isotropic with each decay product having energy m^/S, 
the energy squared of the neutrino is (mg/9)(l + /3^/3), where /? = (1 — mf/m^y/^. 
The branching ration for this decay is 1/9. Thus we obtain 

and 



< Nz^ >t.,^ ( -j (1 + /3V3) = 0.012(1 + /373) (^.22), 



< Nz >t.,^ = 0.037 (^.23), 



which in the hmit /3 — > 1 are close to (A. 8) and (A. 9). 
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Thus in this case a reasonable approximation to the suppression factors is to use 
the Pi with Eq = nit/3. At large values of ms, where such an approach may fail (due 
to the preferential stopping of the higher-energy primary decay neutrinos, such that 
the spectrum is not dominated by these neutrinos I^^l), the t quark final state can 
be neglected compared with the gauge boson final state when calculating event rates. 
Thus we will use (A. 18) - (A. 21) for the b quark final state and the Pi suppression 
factors for the t quark final state. 

SSt ^ h°h° 

RS do not explicitly discuss this case. However, we can easily adapt their results. 
The main decay mode of the Higgs boson when mh < mw (with branching ratio ~ 0.9) 
is to bb pairs. (The Higgs boson final state can in general be neglected compared with 
the gauge boson final states when these are kinematically allowed.) The neutrinos 
occur in the decay of these bb pairs. We can simply regard the decay of the h°h° 
pair as the injection of two bb pairs, with each b quark having a mean energy ms/2. 
This should be a good approximation for ms/2 ^ mb In this case we can use the RS 
results for bb pairs, but with ms —>■ ms/2 and an overall factor of 2. This gives, for 
the non-interacting case, 

< Nz^ >hou^= 3.3 X 10-3(1 + 1.6z^) {A.2A) 

and 

< Nz >ho,.i= 2.8 X 10-2(1 + 0.82z^) (All). 

The suppression factors for the interacting case are [1 — ms/(844GeV)] (z/) and [1 — 
ms/(1.4TeV)] (v) for the < Nz^ >Fi and [1 - ms/(1.5TeV)] (u) and [1 - ms/(2.4TeV)] 
(I7) for the < Nz >Fi. 

We can also use this method to estimate the contribution of the secondary decay 
neutrinos to the < Nz^ >Fi and < Nz >Fi for the case of solar S annihilations to gauge 
boson final states. For the case of S annihilations to a pair of Z bosons, the secondary 
neutrinos come from the decay of the Z's to a bb, cc, or rr pair. (Other lighter quarks 
or leptons are stopped in the Sun prior to their decay and can be neglected '^^l) Thus 
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we can use the results of RS for the case of injection of a pair of b or c quarks or 
T leptons each of energy ms/2. (For the b or c quarks, this will overestimate the 
contribution when ms/2 > or E^, since we are then neglecting the slowing of the 
b and c quarks prior to their decay.) The branching ration for Z decay is 0.15 to a 
b quark pair, 0.12 to a c quark pair, and 0.033 to a r lepton pair '^^l Thus we find, 
using the results of RS, that the contribution of the secondary decays is given by 

< Nz^ >zbz.~ 4.9 X 10-^ (A.26a), 

< Nz^ >zcu^ 1-7 X 10-^ (A266), 

< Nz^ >zru^ 3.8 X 10"^ (yl.26c), 

and 

< Nz >zhu~ 4.2 X 10-^ {A.27a), 

< Nz >zc^« 1.4 X 10-3 (A.276), 

< Nz >zrz.~ 1-7 X 10-=^ (A.27c), 

where, for example, < Nz" >zhu denotes the unsuppressed contribution coming from 
Z decays to b quark pairs. Comparing with the primary Z decays, we find that the 
unsuppressed primary decay contribution to < Nz^ >Fi is about 45 times the secondary 
contribution, and that the unsuppressed primary decay contribution to < Nz >Fi 
is about nine times the secondary contribution. Thus, ignoring suppression of the 
secondary neutrinos, we find from (A. 5) that the primary and secondary decay neutrino 
contributions become comparable at ms ~ 1.4TeV for < Nz^ >zj,, ms ~ 2.2TeV for 
< Nz^ >zv, ms ^ 860TeV for < Nz >z^, and ms 1.3TeV for < Nz >zv 
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Figure Captions 

Figure la S annihilation to h° pairs. 
Figure lb S annihilation to W and Z pairs. 
Figure Ic S annihilation to fermion pairs. 

Figure 2a Thermal relic S scalar density (in units of figh^) as a function of As and 
ms (in units of GeV) for mh equal to 60GeV. 

Figure 2b Thermal rehc S scalar density for mh equal to lOOGeV. 

Figure 2c Thermal relic S scalar density for mh equal to 300GeV. 

Figure 2d Thermal rehc S scalar density for mh equal to 500GeV. 

Figure 3a Ge scattering rate for the case mh = 60GeV. The contours correspond to 

(Js-Ge = a) 1000kg-M-\ b) lOOkg-M-^ , c) 10kg-M-^ d) Ikg-^-^ , e) O.lkg-M-^ 

and f) O.Olkg-M-^ 

Figure 3b Ge scattering rate for the case mh = lOOGcV. 
Figure 3c Ge scattering rate for the case mh = 300GeV. 

Figure 4a Rate of upward-moving muons at the Earth's surface due to neutrinos 
from S annihilation in the Sun, for the case mh = 60GeV. The contours correspond 
to Pdctcctor = a) 2.65xlO-2m-2yr"S b) 10-=^m-2yr"\ c) 10-^m-2yr-\ d) IQ-^m-V"^ 
and e) 10~^m~^yr~^. 

Figure 4b Rate of upward-moving muons at the Earth's surface due to neutrinos from 
S annihilation in the Sun, for the case mh = lOOGeV. 
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Figure 4c Rate of upward-moving muons at the Earth's surface due to neutrinos from 
S annihilation in the Sun, for the case nih = 300GeV. 

Figure 5a Rate of upward-moving muons at the Earth's surface due to neutrinos from 
S annihilation in the Earth, for the case mh = 60GeV. The contours correspond to 
Fdetector = a) 2.65xl0-2m- Vr"\ b) 10-^m-Vr"\ c) 10-^m-^yT-\ d) 10-^m-Vr"^ 
Figure 5b Rate of upward-moving muons at the Earth's surface due to neutrinos from 
S annihilation in the Earth, for the case m^ = lOOGeV. 

Figure 5c Rate of upward-moving muons at the Earth's surface due to neutrinos from 
S annihilation in the Earth, for the case mh = 300GeV. 

Figure 6a Rate of contained events due to neutrinos from S annihilation in the Sun, 
for the case mh = 60GeV. The contours correspond to Fdetector = a) 6.4kt~^yr~^, b) 
lkt-Vi-"\ c) 0.1kt-Vr"\ d) IQ-^kt-Vr"^ and e) IQ-^kf-Vr'^ 
Figure 6b Rate of contained events due to neutrinos from S annihilation in the Sun, 
for the case mh = lOOGeV. 

Figure 7a Rate of contained events due to neutrinos from S annihilation in the Earth, 
for the case mh = GOGeV. The contours correspond to Fdetector = a) 6.4kt~^yr~^, b) 
lkt-Vr"\ c) 0.1kt-Vr"^ d) lO-^kf-V"^ and e) 10-^kt-Vr"^ 
Figure 7b Rate of contained events due to neutrinos from S annihilation in the Earth, 
for the case mh — lOOGeV. 
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Fig. 6(b) 
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Fig. 7(a) 
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Fig. 7(b) 
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